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PREFACE 



This is an interim report on an active research field. It 
contains a comprehensive coverage of those factors which influence 
LSA ignition, provides a useful introduction to the mathematical 
treatment of LSA wave propagation and reaches certain conclusions 
about LSA wave ignition. 

During the period between completion of the manuscript and 
publication of this report, investigation of these problems has 
continued and thus many of the speculations reported here have 
since been confirmed. 
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LASER ABSORPTION WAVE PHENOMENA 

Don E. Harrison, Jr. 

Physics and Chemistry Department 

1.0 INTRODUCTION AND SYMBOLS 

(U) The topics discussed here appear in the literature under a 
variety of headings. Chief attention will be focused upon Laser Supported 
Detonation (LSD) waves and Laser Supported Combustion (LSC) waves. It is 
convenient to group these two phenomena under the more general heading of 
Laser Suppor ted Absorption (LSA) xvaves , which includes both surface init- 
iated and gas initiated LSC and LSD waves as t^ell as plasmons or plasma- 
trons , which describe a related phenomenon, a plasma which does not 
detach from the target surface. The best available summary of current 
thinking on these topics is Volume V of the Proceedings of the October 
1973 LASER EFFECTS/HARDENING CONFERENCE in Monterey, California 1 . 

Warning : Understanding has increased, and many individual positions on 
controversial topics have been modified since its publication! 

(U) The introductory article in reference 1 by P. E. Nielsen 2 is a 
valuable overview of many of the topics discussed here. The entire field 
of study depends to a large extent upon the initial analysis by Raizer 3 , 
published in 1965. A useful discussion of the relationship between the 
various phenomena is contained in a report by Nielsen and Canavan 4 which 
is unpublished. 

(U) This article considers a family of effects which appear when 
laser radiation with irradiance of the order of 10 5 to 10 8 watts/cm 2 is 
incident upon a solid target in air at pressures above approximately 
0.1 atm. Under these conditions a plasma forms at the target surface 
and then may move back along the laser beam if the irradiance is suffic- 
iently large. If the plasma expansion along the beam is subsonic, the 
phenomenon is called an LSC wave; if it is supersonic, the wave is an 
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LSD wave. The general propagation characteristics of the two types of 
waves have been known for some time, but significant aspects of their 
propagation and initiation are still under active investigation. For 
present purposes, initiation has been divided into two steps: ignition 
in which the first few ions and/or electrons are produced, and trans - 
ition which carries the initial ionization into a full blown plasma which 
is propagating as an LSA wave. 

(U) LSA ignition occurs at irradiance levels which are two orders 
of magnitude lower than those normally associated with gas breakdown 
and spark formation, the gas initiated LSA phenomena. The laser-target 
interaction is poorly understood; so a large number of ignition mechanisms 
have been proposed. Hopefully this review can decrease the size of this 
list of candidates. 

(U) Once ignited, the transition to full LSA wave status occurs by 
well understood physical processes. Unfortunately the details of the 
transition depend upon the characteristics of the laser (particularly 
its wavelength) and upon the interaction of a number of competing pro- 
cesses. The role of each process can be identified, but its relative 
contribution to the transition phenomenon is uncertain. As a result, no 
detailed theory of transition for a particular target-laser-atmosphere 
system exists, even though transition times can be estimated with some 
confidence. 

(U) This article surveys the propagation of LSA waves so that the 
differences between LSC and LSD waves can be clearly established. It 
also provides background information about hydrodynamics. 
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1 . 1 SYMBOLS 

The symbols listed here are in the order of their appearance in 
the article. 

X wavelength 

p density subscripts 

m metal 
v vapor 

1 ambient air 

2 shocked air 

v speed of vaporizing surface 

G irradiance 

heat of vaporization of metal 
C specific heat of metal 

T temperature 

j x ion current 

j electron current 

positive ion surface work function 
jz5 electron work function 

I g ionization energy of an atom 

Z s atomic sublimation energy from metal surface 

j vaporization current of neutral atoms from a target surface 

k Boltzmann's constant 

v frequency of radiation 

lifetime of an electron heating cycle 
e electron charge 

N e electron density 
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N 

P 

P 

V 

e 

u 

u* 

h 

T 

D 

q 



j 



c 



p 



c 

a 



v 



c 

K 

e 



critical electron density for opaque plasma 
pressure 

specific volume, that is, volume per unit mass V = 1/p 

specific internal energy 

fluid flow velocity in LAB system 

fluid flow velocity in COM system 

specific enthalpy 

volume element 

shock front propagation velocity 
★ 

- G/P2U2 

★ ★ 

= u 2 P2 = uipx 

Oh/aT) p 

(3e/3T) v 

C /C 
P v 

speed of sound 
absorption coefficient 
thermal conductivity 
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2.0 PLASMA IGNITION 



(U) At least ten ignition mechanisms have been proposed since 1971. 

In some cases the proposal consists of little more than a name, while in 
others a relatively complete theory has been developed. Many of the 
ideas were developed through private communications; so it is exceedingly 
difficult to correctly attribute them to their originators. The situa- 
tion is further complicated by the intermixing of ignition and transition 
processes in the presentation of a particular mechanism. Both Nielsen 2 
and Walters 6 attempt to give proper credit; so these references should 
be consulted by the historian. 

(U) This review will take another direction. The ignition mechanisms 
will be grouped by type, and then the common and distinctive features can 
be more readily identified. The five types are: atom emission, ion emis- 
sion, electron emission, photon mechanisms and shock wave mechanisms. 

This organization will facilitate the description of the ignition pro- 
cess. 

(U) The experimental fact is that a plasma forms near (within 500 
ym) the target surface within a very short time. This time is of the 
order of 50 nsec for LSD waves and somewhat larger for LSC waves. The 
plasma contains ions from both the target and from the atmosphere . In 
all probability it is formed by ions attributable to both sources, but 
there is some evidence that atmospheric ions appear first in the LSD 
environment while target ions appear first in the LSC environment. In 
either event, the transition plasma feeds upon both sources, and the 
plasma which can be examined spectroscopically contains ions from both. 
Because of the nature of this review both possibilities will be discussed 
in some detail; so that the information will be available regardless of 
the ultimate experimental resolution of the present dilemma. The avail- 
able experimental conclusions can then be discussed in context. 

(U) The ignition .of LSA waves is closely related to the more general 
problem of laser-induced particle emission from surfaces. Many experi- 
ments have been performed with ruby lasers (X = 0.6943 ym -*■ 1.78 eV 
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photons) which are quite comparable to the Nd-YAG/glass (X = 1*06 ym 
-►1.17 eV) and C0 2 (X = 10.6 ym -► 0.117 eV) experiments reported by the 
LSA specialists. These investigations have been surveyed by Ready 7 , 
whose book should be consulted for further details . 

Cu) As a general principle, ignition effects rooted in thermal 
processes should vary little between ruby and C0 2 .wavelengths. The 
thermionic emission of electrons must fall into this group. In contrast, 
photon sensitive processes such as the photoelectric effect would hardly 
be comparable because of the difference in photon energy. In the suc- 
ceeding sections ruby laser results will be included when appropriate. 

2.1 IGNITION MECHANISMS 

2.11 (U) Atom emission from the target is one plausible source for the 
plasma if ionization mechanisms are also available. The atoms may be of 
the target species, from a surface oxide layer (in which case the emitted 
species may be a molecule, A10, for example), from surface included 
material (laid down during polishing and finishing) , or from simple 
'dirt* which can find its way to the surface in untold number of ways. 

(U) There will be a characteristic vaporization temperature for 
each of the possible species, a specific heat for each, and a latent 
heat of vaporization for each. Momentum conservation guarantees that 
the emission process would be thermal, which implies that the vaporiza- 
tion temperature determines the velocity distribution of the emitted 
atoms. The mean emission speed of the vapor is determined by the con- 
servation of mass and energy and the parameters of the source and 
emitted material. . 

(u) As an example, consider emission from a perfectly clean metal- 
lic target. Then conservation of mass requires that: 

p v = p v ; 
mm v v 

here v m is the speed with which the vaporization surface penetrates the 
metal, and conservation of energy requires that: 
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v = G/(p [L + C T ]) . 
m ' m 1, in m m J 

The vaporization speed is difficult to determine, so use a relation 
proposed by Anisimov, et.al. 8 : 

v = c exp- (L /C T ). 
m m r v m m nr 

After some algebra the vapor speed, v becomes 

v = AL , (1) 

v m ’ K J 

with 

K 

A = [2£n(p L c /G)] * . 

1 K m mm J1 

The emission speed is of critical importance because it scales the time 
available before the atom leaves the laser beam, and hence the time 
available for the ionization processes. 

(U) In a multispecies environment only a fraction of the irradi- 
ance, G, say ouG, couples to each species, but otherwise the analysis 
is the same, and the emission speed scales with the squareroot of the 
latent heat of vaporization. Actually the problem is even more compli- 
cated because surface species will be heated first; so the steady state 
emission process envisioned here is an idealization. 

(u) Ignition mechanisms of this type must be strongly target 
material sensitive. They will also be sensitive to the target surface 
condition; so the ignition characteristics might change from pulse to 
pulse as the laser flux cleaned the surface. This pulse to pulse vari- 
ation is plausible whether the actual mechanism depend upon an oxide 
coating, inclusions, or dirt. If the laser pulse rate is not too high 
an oxidized surface maintained by the atmosphere, again A1 is a good 
example, would have time to reconstitute between pulses. 

(U) Because this mechanism depends upon pure heating it would 
depend strongly upon the laser intensity, but rather weakly upon the 
laser frequency. 
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(U) Very high energy neutral particles are also emitted from 
targets 7 . The energy of CO molecules have been measured by time-of- 
f light techniques to be 14 eV at 10 MW/cm 2 irradiance and 500 eV at 
140 MW/cm 2 from Q-switched ruby lasers. Thermal emission at these 
energies would require a surface temperature of the order of 10 5 K, 
which is highly unlikely, because the CO molecule would certainly dis- 
sociate at such temperatures. It is possible that the neutral was formed 
by recombination outside the target after emission, or that these high 
energy neutrals were formed by charge exchange. These effects have been 
included here for completeness, but there is only a low probability that 
these high energy neutral species are directly emitted from the target; 
so they are not candidates for ignition source. 

2.12 (u) Ion emission adds an additional complication to the atom emis- 
sion discussed above. If the surface temperature is sufficiently high 
an appreciable fraction of the emitted atoms will be ionized. If the 
emitted particles are singly ionized the Richardson-Smith equation 
applies, and the' ion curre'nt is given by: 

j + = A + T 2 exp(-0 + /kT) 

where A is a constant, T is the absolute temperature, and 0 is the 
positive ion work function. 0 + is essentially unknown, but it can be 
estimated by the argument that the energy required to remove the singly 
charged ion, 0 + , plus the energy required to remove a single electron, 
should equal the energy required to remove a neutral atom, plus 
the ionization energy of the neutral atom, I . This argument applied 
to A1 gives 0 + (Al) = 5.3 eV, and to Cu gives 0 + (Cu) = 7.25 eV. For 
comparison, 0(A1) = 3.6 eV and 0(Cu) = 4.4 eV. Thus the positive ion 
emission will be considerably smaller than the electron emission at the 
same temperature. The Richardson-Smith equation is only an approxima- 
tion; so it should be used with caution. 

(U) Even less reliable is the Saha-Langmuir equation which is 
derived under equilibrium assumption, but is often applied to determine 
the fraction of the particle current which is ionized by the laser beam. 
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In particular, target surface temperatures obtained in this way are 
highly suspect. The Saha -Langmuir equation predicts that the fraction 
ionized takes the form 

j + /j 0 ~exp[(tf-I a )/kT]. 

As the ionization potential is generally larger than the work function 
the degree of ionization must increase rapidly with temperature. From 
Q-switched ruby laser data at 6 MW/cm 2 the duration of the emitted ion 
pulse is ~ 100 nsec, which is of the order of magnitude required for 
the prompt ignition of LSD waves. 

(U) Time-of-f light energy measurements on the ions emitted by 
Q-switched ruby laser experiments give energies of the order of 1000 eV. 
As surface temperatures never approach such values plasma effects must 
be influencing the ion emission process. Multiphoton effects of this 
magnitude are considered to be totally out of the question. 

(U) Ion emission from oxide layers and/or 'dirt' is certainly pos- 
sible, but is too dependent upon the details of a particular targets 
condition for general discussion. High energy molecules containing atoms 
from surface species, particularly oxides, are readily detected as 
molecular ions, but the ionization may have occurred after emission. 

(u) The preceding discussion makes use of temperatures, which im- 
plies that the materials are at least approximately in thermal equilib- 
rium. This assumption may not be justified, even though no non-equilib- 
rium effects have ever been positively identified. In particular, in 
the presence of the electron emission processes discussed in the next 
section, strong electrostatic forces must be applied to ions near the 
target surface, and non-equilibrium emission processes become possible. 
When many electrons have been stripped from the surface region, heavy 
positively charged ions remain; so that a large fraction of the target 
atoms emitted could be ionized. Charge conservation at the target sur- 
face becomes a major consideration. 

(U) In summary, it is possible to visualize processes by which the 
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target could emit a significant number of ions in a short period of time, 
but the processes depend upon detailed considerations of particular mater 
ials and surface conditions which defy analysis at the present time. 

2.13 (U) Electron emission processes are the most likely ignition mech- 
anism for LSD waves. Thermionic emission is known to be present, and if 
the laser pulse shape is known the shape of the electron pulse can be 
calculated from the Richardson-Dushman equation. This equation can be 
derived from either thermodynamic 9 or kinetic 10 theory arguments. The 
usual form is 

j_ = A_T 2 exp(-ji/kT) t 

where j is the electron density, <t> the electron work function of the 
target surface, and A is a universal constant (A = 60.2 A cm 2 deg 2 ) 
for many metals. Ready 7 shows that the Richardson-Dushman equation is 
a reasonable description of the electron emission from normal laser 
pulses, and that thermionic emission of the order of mA/cm 2 is readily 
produced within 50 nsec by Q-switched ruby lasers 11 . These measurements 
were performed at 10~ 8 torr. At higher pressures where absorbed gases 
become more prevalent the electron emission rises, which indicates that 
electron emission can also occur from surface materials. 

(U) Thermionic emission is also influenced by surface morphology. 

If the surface contains pits, protrusions, or inclusions, the local 
temperature may be much higher than would be predicted from the laser 
flux and thermal properties of the bulk target material. The tempera- 
ture at the tip of a protrusion will be significantly higher than that 
of the surrounding material because the tip is unable to dissipate heat 
to the bulk of the material so rapidly as a planar surface can. 

(u) Surface atoms , whether oxide layers, or 'dirt', disrupt the 
lattice periodicity of the target, and affect the work function of the 
surface. At 0 = 5 eV, a one percent reduction of the work function 
causes an 80 percent increase in j at 1000 °K. A few surface atoms 
produce noticable effects, a monolayer of coverage completely changes 
the work function. Because all materials emit thermionically , at high 
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irradiance, target material effects will appear primarily in the time 
dependence of j , that is in the shape of the electron pulse ignited by 
the laser pulse. 

(U) ' Field enhanced emission from cold targets is a very prompt 
mechanism, and therefore is attractive as a LSD wave ignition process. 

Musal has explored this possibility 12 from the viewpoint of the Fowler 
and Nordheim analysis 13 which is also discussed in the Dushman article 
cited above 9 . Musal has done the calculations for laser irradiances in 
the range from SO to 500 MW/cm 2 and finds that the electric field 
strength is approximately two orders of magnitude too small for the effect 
to be significant for electron emission from a plane surface. However, 
when surface features are included in the analysis the field strength is 
greatly enhanced, and the mechanism becomes feasible. Photon assisted 
field emission has recently been discussed by Caroli, et.al. 14 . 

(U) Schottky emission , thermionic emission in the presence of an 
external electric field, has also been extensively studied. This is a 
quantum mechanical barrier penetration process in which the work function 
is field dependent 15 . The analysis shows that the work function will be 
reduced, but that otherwise the process is similar to thermionic emission 16 
Fields of the order of 10 5 volts/cm cause a very small lowering of d. The 
power of T in the pre-exponential factor of the Richardson-Dushman equa- 
tion is difficult to determine experimentally, particularly in the field- 
enhanced environment. The factor may well be T 3 rather than Richardson's 
T 2 in the LSA experiments, because there is considerable experimental 
simularity between the ignition experiments and the 'complete photo- 
electric effect' 16 . 

(U) Multiphoton emission processes can produce significant electron 
current density from ruby and Nd doped glass laser photons. The mechanism 
is hardly feasible from the low energy photons of the CO 2 laser. The 
mechanism has been included here because it gives the earliest response 
of any of the mechanisms proposed. Multiphoton processes become partic- 
ularly significant for pico-second laser pulses when the irradiance is 
small so that there is little temperature rise in the target. Thermionic 
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emission decreases exponentially as the irradiance decreases, while a 
two-photon emission process decreases as the square of the irradiance; 
so at some level the two effects must be equal, but this normally occurs 
at irradiance levels well below those required to initiate LSA waves. 

2.14 (U) Photon emission mechanisms are poor candidates for LSA ignition 
The two immediate possibilities are reflected photons from the laser 
pulse and photons generated in the target. The reflected photons are not 
sufficiently energetic to cause appreciable ionization, but their effect 
can be included under the field emission heading because they do increase 
the effective electric field 12 . Photons from the target could be either 
thermally emitted by the hot target surface, or emitted as the result of 
a multiphoton atomic process which occurs at the surface. 

(U) Although sufficient energy is available to melt the target sur- 
face and/or vaporize it under some conditions, the temperatui'es reached 
are realtively low, and only low energy photons would be produced. Such 
photons probably exist, but they can have little effect upon the transi- 
tion from ignition to plasma; so it seems unlikely that they can be major 
contributors to the LSA initiation process. 

(U) Multiphoton atomic process could make some contribution for 
ruby laser or Nd/glass illumination, but not for CO 2 illumination because 
the photon energy is so small that 10-20 photon process would be re- 
quired. There is not experimental evidence that photon processes contrib 
ute significantly to LSA ignition. 

2.15 QJ) Thermomechanical response of the target to the laser pulse has 
been proposed as an ignition mechanism by. Steverding 12 . In this model, 
rapid heating in a shallow region near the target surface produces a 
compression stress wave which propagates to the target surface as a 
shock wave. The shock wave produces spallation at the surface. The 
emission velocity is the shock velocity, and the emitted material is a 
plasma because the shock wave has both vaporized and ionized the target 
material. 
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(U) As envisioned, the target is heated for a relatively long time 
(20 ysec) and the laser energy is deposited in a very shallow region 
(10 -1+ cm = 1 ym) . This last assumption requires a depth of penetration, 
a skin depth, of approximately X/10 for the CO 2 laser pulse. A region 
of this depth is consistent with the calculation of Ready 1 *, for the 
first 50 nsec of the pulse, but thermal conduction increases the depth 
by an order of magnitude within the first 200 nsec. This mechanism 
appears to be too slow for LSD ignition, but it may be feasible for LSC 
ignition. 

2.2 IGNITION EXPERIMENTS 

(U) Experimental study of LSAwave ignition processes cannot pres- 
ently be separated from effects attributable to transition processes, 
but certain mechanisms can be ruled out by more general properties of 
LSA waves. Here some of the experimental properties of LSC and LSD 
waves are summarized and interpreted in the light of the various pro- 
posed mechanisms. Ignition of LSC waves differs sufficiently from LSD 
ignition that they are considered separately. Much more experimental 
effort has been devoted to LSD waves, see the reports by Hall, et.al. 18 ; 
Walters 6 ; and Musal 12 ; so the ignition properties are more completely 
determined. Pertinent LSC results have been reported by Fowler and Smith 19 
Conrad, Mangum, and Gurley 20 ; and Thomas 21 who has used computer simula- 
tions to good effect. 

(U) This survey will not discuss LSB (Laser Supported Blast ) 
effects in detail. These spherical shock waves occur in conjunction with 
LSC and LSD waves, but are not supported by the laser. Hall, et.al. 18 
have pictures and a discussion. Target vapor 'microjets' are also 
observed 18 , which move normal to the target surface; so if the target 
is inclined with respect to the laser beam direction they are easily 
distinguished from LSC and LSD waves. Micro jets ’ignite* much more 
slowly than LSD waves, but on a time scale comparable to LSC wave 
ignition times. 
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(U) Vacuum ignition of A1 targets by nsec duration Nd-glass laser 
pulses of irradiance ^ 10 19 W/cm 2 have been studied in some detail 
by Brooks 22 . Note that the targets were held in a relatively good 
vacuum 10 3 torr) , and the irradiance is three orders of magnitude 
greater than that associated with LSD wave ignition, but that the total 
energy absorbed by the target is comparable to that reported in the LSC 
and LSD experiments . 

2.21 (U) Laser supported combustion wave ignition occurs in the irradi- 
ance range 2 x 10 3 < £ 6 x 10 9 W/cm 2 . Hall, et.al. 19 has observed 

LSC waves in the short pulse (~ 20 ysec) CO 2 laser system in which their 
LSD experiments were performed. They have not made extensive studies 
of LSC wave initiation or propagation. The detailed information currently 
available comes from the cw CO 2 laser studies of Fowler and Smith 19 , and 
the long pulse 1 sec) CO 2 laser studies of Conrad, Mangum, and Gurley 29 . 

(U) When compared with the LSD wave experiments the single most 

outstanding difference is the relatively long initiation time of the LSC 

waves. If the initiation time is defined as the time until the LSA 

becomes detectable, both ignition and transition must occur during the 

initiation period. This time is of the order of 200 ysec for LSC waves, 

which must be compared with LSD initiation times of the order of 50 nsec. 

LSC is very much slower, more than three orders of magnitude, so the 

ignition processes must be very different. On the assumption of a linear 

rise and that G occurs after 200 ysec, the target surface absorbs 
P 

~ 100 J/cm 2 from a G = 10 6 W/cm 2 laser before initiation. With this 

P 

time span and this much energy available, even the slower ignition mech- 
anism, such as atom and ion emission, are feasible. When comparable 
laser irradiations are applied to targets in vacuum chambers, where the 
pressure is of the order of 0.1 atm, plasma micro jet phenomena are 
observed. Thus target vaporization with consequent emission of atoms 
and/or ions seems to be an observable characteristic of this irradiance 
region. Thomas 21 has- done computer simulations based upon vaporization 
ignition which compare very favorable with the MSNW experiments 23 . 
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CU) In summary, target vaporization mechanisms appear to be suffic- 
iently to explain the existing information about LSC ignition, but the 
necessity of these mechanisms has not been established. The faster 
electron emission mechanisms certainly have the opportunity to function, 
but the irradiance may be insufficient to initiate a surface plasma. 
Sufficient energy has been dumped into the surface that thermionic 
emission must occur but the field strength contribution from the irradi- 
ance minimized the Schottky emission contribution. Thus there is indirect 
evidence that the Schottky mechanism is an important contributor to LSD 
ignition. 



2.22 (U) Laser supported detonation wave ignition occurs in the irradi- 
ance range 2 x 10 7 <_ <_ 5 x 10 7 W/cm 2 . Most of the pertinent experi- 

mental information on LSD ignition has been obtained by Hall, et.al. 18 , 
and by Walters 6 who has also done a detailed integration of the experi- 
mental evidence applicable to the ignition problem available at the time 
of his review. 



Cu) LSD waves ignite very promptly; the initation time is of the 
order of 50 nsec. Quite independent of the target material the initiation 
occurs after the target has absorbed approximately 1.7 J/cm 2 of energy. 
This value was obtained by Walters 6 who integrated his experimental laser 
pulse shape out to the initiation time. This is insufficient energy for 
uniform surface vaporization to occur, but it is consistent with either 
a thermionic or Schottky emission mechanism. At this irradiance the rms 
electric field strength is insufficient to produce field emission from a 
plane surface; so the Schottky mechanism is indicated. 

(U) Walter's careful studies of the ignition process indicate that 
LSD waves ignite from 'numerous individual luminosity sites'. The wave 
front forms when the emission from these sites combine in front of the 
target surface. He has done careful scanning electron microscope (SEM) 
studies of target surface, and has compared the appearance of an indi- 
vidual surface feature before and after LSD wave ignition. The pictures 
are in ref. 6. A feature identified before the laser pulse is immediately 
recognizable after the pulse. The SEM scans an area 50 by 250 ym at a 
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2000X magnification. Although absolute correlation between surface 
features and luminosity sites has not been accomplished the evidence 
that features cause luminosity sites is very strong. 

(U) The laser pulse does cause surface damage, and it does change 
the surface features, but there is no gross change in any feature which 
could be associated with significant vaporization of its surface. When 
all of the data are taken together electron emission from surface fea- 
tures appears to be the most likely ignition mechanism. The ignition 
may be thermionic, field emission as proposed and discussed by Musal* 2 , 
or a Schottky effect combination of the two at the surface feature. 
Vaporization from surface features can certainly occur, but it must be 
a second order effect because the features exhibit so little damage, 
and electron emission is dominant. 

(U) Tables I and II from Hall, et.al. 16 give LSD ignition thresholds 
as measured by image convertor photographs and by measurement of a re- 
flected signal. The agreement is relatively good when the lack of 
uniformity of the surfaces is considered. 

(U) Table III is Walter's summary 6 of the relative importance of 
the various LSD ignition models which have been proposed. 
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TABLE I 

THRESHOLDS FOR IGNITION OF LASER-SUPPORTED ABSORPTION WAVES USING 
A 28 CM FOCAL LENGTH LENS AND USING A NEW SPECIMEN FOR EACH TEST 
AT ONE ATMOSPHERE. THRESHOLD DETERMINED BY IMAGE CONVERTER PHOTO- 
GRAPHS OF LUMINOSITY. 



Material 


E,p(j oules/cm 2 ) 


%iax^ watts ^ CI1 


Copper Foil 


266 


3.83 


X 


10? 


Lucite 


258 


3.72 


X 


10 7 


Polished Tungsten 


229 


3.30 


X 


10 7 


Copper Plate 


192 


2.8 


X 


10 7 


Tungsten Plate 


186 


2.68 


X 


10 7 


(unpolished) 


Aluminum 


176 


2.53 


X 


10 7 


(7075 alloy) 


Aluminum Foil 


160 


2.30 


X 


10 7 


Tantalum 


128 


1.84 


X 


10 7 


Carbon 


88 


1.27 


X 


10 7 


Aluminum Black 


87 


1.25 


X 


10 7 


(very heavily anodized Al) 


Lead Foil 


64 


.92 


X 


10 7 


(polished) 


Lead Plate 


43 


.62 


X 


10 7 


(unpolished) 
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TABLE II 

THRESHOLDS FOR IGNITION OF LASER-SUPPORTED ABSORPTION WAVES USING 
A 28 CM FOCAL LENGTH LENS AND USING A NEW SPECIMEN FOR EACH TEST 
AT ONE ATMOSPHERE THRESHOLD DETERMINED BY MEASURING REFLECTED 
SIGNAL. 



Material 
Fused Silica 
Lucite 

Lexan (polycarbonate) 

h 2 o 

Teflon 

Copper 

Painted Aluminum 

(.001 cm gloss black lacquer) 

Titanium 

Fiberglass Epoxy 

Cork 

Nickel 

Tungsten 

Soda Glass 

Painted Aluminum 
(.004 cm red enamel) 

Aluminum Foil 

Aluminum 7075 Plate 

Lead 

Alumina 



T (joules/cm 2 ) q (watts/cm 2 ) 
i max 



310 


4.5 X 10 7 


310 


4.5 X 10 7 


280 


4.0 X 10 7 


220 


3.2 X 10 7 


195 


2.8 X 10 7 


191 


2.75X 10 7 


173 


2.5 X 10 7 


160 


2.3 X 10 7 


158 


2.3 X 10 7 


144 


2.1 X 10 7 


137 


1.97 X 10 7 


124 


1.79 X 10 7 


106 


1.52 X 10 7 


85 


1.22 X 10 7 


84 


1.2 X 10 7 


75 


1.08 X 10 7 


68 


.98 X 10 7 


59 


.85 X 10 7 


42 


.60 X 10 7 



Painted Aluminum 

(.004 cm flat black lacquer) 
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TABLE III 

ASSESSMENT OF THE ROLE OF VARIOUS MECHANISMS IN LSD WAVE 
INITIATION ON PRACTICAL ALUMINUM SURFACES 





Mechanism 


Role 


Experimental Basis 



(1) Uniform target 
vaporization 


None 


SEM, photography, breakdown 
time 


(2) Oxide vapor 
absorption 


Minor 


Oxide-free surface initiation 


(3) Chemical reaction 


Unknown* 




(4) Thermionic emission 


Ma j or 


Target emission, breakdown 
time, pressure dependence 
of threshold 


(5) Planar reflection 


No inde- 
pendent 
role for 

G <5X10 8 
P 


Clean-air breakdown 


(6) Local target heating 


Major 


SEM, photography 


(7) Contaminant and gas 
desorption 


Variable 


Contaminant -free surface 


(8) Field emission 


Minor for 
GP<3X10 8 


SEM, bent sample, breakdown 
time 


(9) Shock heating 
of the air 


Unknown* 




( 1 0) Nonequi 1 ibr ium 
ionization of 
vapor 


Unknown* 




(11) Defect -enhanced 
fields 


Minor for 

G <3X10 8 
P 


Breakdown time 



* Spectroscopy now underway should provide assessment. 
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2.23 (u) Plasma emission experiments are discussed in this section. 
Although there is no direct evidence that these results bear directly 
upon the LSA ignition problem the similarities are sufficiently striking 
that the connection may ultimately be made. 

(U) The experiments, performed by Brooks 22 , used a Q switched Nd- 
glass laser with energy between 2.5 and 12 joules. The pulse duration 
was 25 nsec, and the focal spot size was '*'0.018 cm 2 ; so the irradiance 
was 1.4 x 10 9 <_ Gp 5.7 x 10 9 W/cm 2 . Note that the total energy fluence 
absorbed by the target (140-700 J/cm 2 ) is comparable to that reported in 
the LSC and LSD experiments . 

(u) The laser pulse produces a burst of plasma from an A1 target. 
Careful analysis of the plasma by probes determined that it consists of 
3 components which are described in Table IV. 

TABLE IV 

Emission 

Component Time Speed Density 

first early 3 nsec 1.1 x 10 6 m/sec 2.5 x 10 12 /cm 3 

second early 12 nsec 5.9 x 10 5 m/sec 2.5 x 10 12 /cm 3 

main 0-10 nsec 1.1 x 10 5 m/sec lO^/cm 3 

(U) A very prompt (< 1 nsec) electromagnetic pulse was also detec- 
ted. Experimental studies confirmed that its wavelength is less than 
0.4 y; so it was not a reflection of the incident laser pulse. 

(U) An applied magnetic field has little effect on the early pulses; 
so they cannot be purely electronic in nature. Heavy ions must be present 
hence the term "plasma pulse" seems appropriate. The main pulse is delaye 
by the applied magnetic field (B = 950 gauss) . The delay is consistent 
with an interpretation which finds the plasma front relatively rich in 
electrons while the main body is richer in ions 21 * 26 . 

(U) The early plasma pulse (s) have been mapped with probes. Figure 
1 shows the plasma extent after 20 nsec. The two early pulses have not 
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AXIAL POSITION (CM) 



Fig. 1 Rftlativfi density contours for a laser produced 
plasma 20 nsec after ignition. :!ee text for details 



20A 



UNCLASSIFIED 



separated at this time. Note that the plasma is asymmetrical and appears 
to be expanding along the specular reflection direction. The plasma has 
expanded out to the beam spot, and is now eroding the entire target sur- 
face. Experiments performed on large area targets show that after 50 
shots the naked eye visual erosion region has an area of 1 cm 2 when the 
focal area is 0.018 cm 2 . Microscopic analysis of a larger region has not 
been carried out, Figure 1 suggests that the plasma extends beyond this 
range . 

(u) The plasma contours after 160 nsec are shown in Figure 2. 

Note that the two early pulses have now separated from the main plasma 
(one is still visible on the right hand side of the figure). The plasma 
emission is now symmetrical, and the main pulse is moving normal to the 
target surface. 

(u) The situation after 800 nsec is shown in Figure 3. Note in 
particular that the plasma is still attached to the target; so that it 
is still feeding upon the target by sputtering. 

(u) This material has been included because it demonstrates the 
target material behavior when it expands into an atmosphere of comparable 
density (n CT 10 13 /cm 3 ). At a standard atmosphere the plasma would inter 
act much more strongly with the background gas; so the development must 
be much more complicated than illustrated here. Still certain similari-' 
ties should remain. The prompt emission from the target should not be 
background pressure sensitive. The electron-rich front surface of the 
main plasma burst is a target phenomenon which should carry over. Note 
that it supports the idea that the initial emission must be electronic 
with ions following closely. This effect will be discussed in greater 
detail below. Even though the laser pulse has been extinguished for 
almost 0.8 ysec in Figure 3 the plasma is still feeding on the target 
as it expands away from the surface. There is strong evidence that the 
radial expansion of the plasma is constrained by a self-induced magnetic 
field 22 " 27 . This would help to explain the long time of contact between 
the plasma and the target surface. 
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AXIAL POSITION (CM) 



Fig. 2 Relative density contours for a laser produced 
plasma 160 nsec after ignition. 
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(U) Electron microscope examination of the focal area on the tar- 
get produces evidence for shock wave damage of the target. "Large" 
masses of A1 have been spalled from the target. There is also evidence 
for "splashing" of liquid A1 from the target onto the probes when they 
are located close to the target surface. At these very high irradiances 
and with these very short pulses, Steverding's model 17 seems to apply 
during some part of the ignition event. 

(y) These investigations also have a bearing upon the nature of 
the electron emission mechanism which causes prompt ignition. The 
electron energies corresponding to the emission speeds of the plasma 
components, see Table IV, are first early =3.4 eV, second early = 1 eV, 
main = 3 x 10 2 eV. These results are consistent with an interpretation 
which suggests that the first two components are caused by a high field 
emission mechanism, while the main component is thermionic. This inter- 
pretation has been discussed by Thomas, Musal and Chou 28 . The structure 
of the components suggests that the electrons are emitted first and that 
ions from the target surface are then pulled along behind the electrons 
by electrostatic forces. 
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3.0 QUALITATIVE SUMMARY OF TRANSITION PROCESSES 



(U) The ignition process is complete as soon as conditions have been 
established in the air in front of the target surface which will ultimately 
lead to the production of an LSA wave. A certain number of charged par- 
ticles must be present (both electrons and ions) and the number of both 
must increase with time. The transition to full LSA behavior is extremely 
complex because a large number of processes, both productive and dissi- 
pative, are occurring simultaneously. If an LSA wave is to result, the 
productive mechanisms must dominate. Any mechanism which increases the 
absolute number or density of electrons is productive, any which decrea- 
ses the electron population is dissipative. Attention can be focused on 
the electrons because an electron produced is also an ion produced. The 
converse of this statement is not necessarily true because electrons can 
diffuse out of the ionized region when it is small. In most cases this 
surface leakage is not a dominant consideration, but it is a significant 
contributor when the dimensions of the ionized volume approach the mean 
free path of electrons in air. This number depends upon many considera- 
tions, but is of the order of 0.5 ym in STP air. In general focal spots, 
and hence ignition plasmas are much larger; so electron diffusion is a 
secondary consideration. 

(u) The atomic-scale production and dissipation processes are sum- 
marized in Figure 4. In the following material A will stand for an 
unexcited atom or molecule. A* for an atom in an excited state and A+ 

for the singly charged ion of that atom. Multiple ions, A 2 + etc., and 

the transition from the molecular to the atomic state will be ignored in 
this discussion. Both complicate the already complex discussion without 
significantly contributing to the analysis. The atom-molecule transition 
is an additional loss mechanism which reduces the effectiveness of the 
production mechanisms discussed below. 

(U) All of the productive processes feed directly or indirectly upon 
the energy which the laser pulse supplies to the ionized region. This 

region will be referred to as the plasma, even though in the strictest 
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X + ( ) = X ' 

WHERE X = A, A* OR A + 



( _ e - ) 




( + e-) 



Fie. 4 Schematic representation of production and 
dissipation processes. Arrows to the right represent 
an ir>crfip?A in ionization. Arrows to the left rep- 
resent a decrease in ionization. 
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sense absolute charge neutrality can not be guaranteed because of electron 
diffusion and target surface effects. 

(U) Because photons couple poorly to atomic mass particles (the mass 
mismatch is too large) photon-electron processes play the primary role. 
These processes can be catagorized as bound-bound, bound-free, and free- 
free where the descriptive words apply to the condition of the electron 
before, and then after the photon- electron collision. Although these 
processes must occur 'near' an atom or ion to conserve momentum, the 
photon-electron collision mechanism dominates the event. The processes 
are 

1) hv + A -»• A* (bound-bound) 

2) hv + A (A*) A + + e (bound-free) 

_ I 

3) hv + e -> e (free-free) 

Mechanisms 1) and 2) combine to produce multiphoton ionization , or on a 

macroscopic scale, microwave heating of the plasma 28 . The two methods of 
analysis are formally equivalent. This is a plausible plasma heating 
mechanism where the laser photon energy exceeds 1.0 eV (Ruby, Nd-Glass, 
Nd-YAG) but it is not practical for the C0 2 laser photon, which is 
approximately 0.1 eV. 

(U) In the absence of multiphoton ionization, mechanism 3) must be 
the dominant coupling between the laser pulse and the plasma. 

(U) Given a supply of energetic electron from mechanism 3) , addit- 
ional ionization can be produced by electron-atom mechanisms. The 
productive mechanism is: 

4) e - + A(A*) -> A + + 2e . (ionization) 

Note that ionization will be facilitated by 

_ r 

5) e~ + A -> A* + e . (excitation) 

Mechanism 4) is the ultimate plasma source, while in one sense mechanism 
5) is a loss mechanism because it withdraws energy from the electron 
heating process. It's exact catagorization will depend upon the details 
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of a particular plasma condition and upon the electron and ion tempera- 
tures in that plasma. 

(if) Two additional mechanisms: 

_ _ f 

6) e + A* -*■ A + e (collision de-excitation) 

and 

7) e + A A + hv' (recombination) 

are clearly dissipative. 

(U) A great deal of study has been done on photon-electron and 
electron-heavy particle scattering. Much of the effort required to trans- 
late these investigations into forms useful for the LSA transition analysis 
has been done by Ready 7 (see in particular the sections on spark ignition 
in air) and by Zel 'dovich and Raizer 2 9 (particularly in Vol. I, which 
contains the fundamental physical material). 

(If) It should be noted that a description of the transition from the 
atomic interaction mechanisms to plasma ignition is dependent upon some 
form of transport theory, which in turn requires local (psuedo) thermo- 
dynamic equilibrium. The electron equilibration time is such that a well 
defined electron temperature seems plausible, but the concept of an ion 
temperature should be used with great caution, because much of the ions' 
energy may be associated with the streaming of the plasma away from the 
target. Also, the ion equilibration time is approximately four orders 
of magnitude greater than the electron equilibration time, which places 
it uncomfortably close to the entire initiation time of the LSA wave. 

In summary, the electrons probably have a temperature, the ions may have 
a temperature, but the electrons and ions are clearly not in thermal 
equilibirum. 

(U) To a good approximation the plasma density will grow exponentially 
with the number of high energy electron "generations" n, required to 
reach a given level of ionization. Thus 
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where T c is the lifetime of one electron heating cycle: 

_ i 

hv + e -*■ e , 



followed by 

€ + A (A*) * A + + e' . 

For an ignition plasma of the type studied by Brooks 22 approximately 16 
generations are required to produce plasma densities approaching 100% 
ionization. If N q were 10 electrons, 43 generations would be required. 
99% of the ionization occurs in the last 7 generations. The breakdown 
process can proceed in times as short as a few nseconds 29 . 
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4.0 PROPAGATION OF LSA WAVES 



(u) The preceding two sections discuss the ignition of a plasma at 
the target surface and the transition of that plasma to a high density 
state. This section is devoted to the behavior of the system once the 
plasma has become dense. 

(U) The transition process is complete once the plasma reaches a 
constant density. If this density is less than that established by setting 
the plasma frequency equal to the laser frequency 

N e 2 

e 

= v = V- 

m p laser 

v = 8.97 x 10 3 N *” (cm -3 ) 



or 

N p (cm -3 ) = 1.12 x 10 9 /X 2 (m) 

the plasma is transparent to the radiation, and the target surface continues 
to supply material to the plasma. Once the plasma density exceeds the 
threshold value the plasma becomes opaque and shields the target from the 
laser pulse. 

(U) Three types of behavior have been identified experimentally. When 
the irradiance, G, is of the order of 10 4 W/cm 2 a plasma is formed at the 
target surface, but its density never exceeds the plasma frequency threshold. 

The plasma is in a steady state under the combined influences of the laser 
irradiation, convection, reradiation and target surface erosion. The 
resultant PLASMATRON is transparent to the laser radiation; so it remains 
attached to the target surface. When the irradiance is of the order of 
10 5 W/cm 2 the plasma detaches from the target surface and propagates rather 
slowly (v - 10 3 cm/s) back up the laser pulse. The plasma is opaque and 
protects the target surface from the laser beam. The ignition mechanism 
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in the air in front of the plasma is very complex. This low velocity dis- 
turbance is the LSC wave. Lastly, when the irradiance is of the order of 
10 7 W/cm 2 , the plasma propagates with a velocity which exceeds the speed 
of sound and an LSD wave results. The names of the waves come from their 
similarity to high temperature hydrodynamic waves. The same mathematical 
models apply; so hydrodynamic theory can be used to describe LSC and LSD 
waves. This analysis was first developed in the Soviet Union by Raizer 3 
and in the West by Ramsden and Savic 3 ®. 

(U) Detonation waves are easier to discuss theoretically because a 
number of simplifying assumptions can be made in the analysis. Combustion 
waves are more complex and the details of their propagation are less well 
understood. The analysis here passes from a simple one-dimensional dis- 
cussion of the relevent hydrodynamic equations to a presentation of the 
salient characteristics of the LSD wave, and then considers the LSC wave 
where closed form analytic expressions are less useful. 

4.1 HYDRODYNAMIC EQUATIONS 

(U) Consider a one-dimensional system of the form shown in Figure 5, 
which is characteristic of the LSD wave. Here the absorption region is 
optically thick; so the laser radiation is absorbed there. The wave prop- 
agates to the right with speed D. The region denoted "ambient air", will 
be Region 1 in what follows, and is "ahead of" the wave. The region 
denoted "relaxation region" will be Region 2, and is "behind" the wave. 

(U) The analysis presented here is not sophisticated, but it includes 
all of the physical processes and variables which are pertinent to this 
discussion. It differs from the similar hydrodynamic derivations because 
an energy source term has been included, but the momentum contribution 
from the laser pulse have been neglected. A more detailed discussion is 
available in Zel'dovich and Raizer 31 . In the material which follows the 
thickness of the discontinuity (•£ in Figure 5) will be neglected. 

(U) The variables required to characterize the system are: 
p the p ressure 
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Fig. 5 Temperature and pressure profiles typical of a 
laser- supported-detonation wave. 
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p the density or its reciprocal 

V = 1/p the specific volume 
e the specific internal energy 

u the fluid flow velocity in the LAB system 

Occasionally, h = e + PV the specific enthalpy will be used to simplify 
equations. Here, the term "specific" means "per unit mass". The general 
term "fluid" has been used rather than "air" or "gas" because the material 
behind the discontinuity is a plasma. 

(U) The analysis considers the conservation of an element of mass. 

Am, its momentum, and its energy, as the energy, as the element moves 
across the discontinuity. In terms of the variables, define a volume 
element, St, on each side of the discontinuity. Then, 

St = (D - ui) • ASt = (D - U 2 ) * ASt; 

where the area of the element A, remains constant as it passes through the 
discontinuity, and of course the same time duration, St, applies. 

(U) Then 

Am = pSr 

Am = p 1 (D - ui) ♦ ASt = P 2 (D - U 2 ) • ASt; 

which gives 

Pl (D - ui) = p 2 CD - u 2 ) • 

This conservation of mass relation is most conveniently written as 

D(P2 - Pi) = P2^2 ‘ Piui* I' 

(u) Conservation of linear momentum requires that the momentum change 
of the mass element be equal to the impulse received. Then 

(P 2 - Pi)A5t = Amu 2 - Amu!. 

•A - -V ^ _ \ 

When D, A, ui, and U 2 are all parallel this equation simplifies to 

D(p 2 u 2 - Piui) = (P 2 + P 2 U 2 ) - (Pi + Piu 2 ). II' 

2 1 

This is the simplification which reduces the relation to a one-dimensional 
form. 
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(U) Conservation of energy follows from the First Law of thermo- 
dynamics : 



AU = AQ + AW. 



The element of mass, and hence the quantity, U, contains energy in two 
forms: the specific internal energy , e, and the specific kinetic energy 
(1/2) u 2 . (Watch the units!). Then one writes 

AU = Am (e2 - ei) 



where 



e = e + (l/2)u z . 



The heat added is the integrated laser flux through A in time fit, 

AQ = G Afit 
o 

and the work done on the element of mass is 
AW = (P 2 U 2 - Piui) » Afit. 



When all terms are combined the relation reduces to 



D(p 2 £ 2 -Piei) = G + (p 2 u 2 £2 + P 2 U 2 ) - (Piui£i + Piui). Ill' 

(U) These equations can be considerably simplified by transforming 
to a coordinate system which moves with the discontinuity: 

u* = u - D. 



Note that e / e* because of the kinetic energy term. The three conser- 
vation equations reduce to: 

Piu* = p 2 u* s j, I 

Pi + Piu* 2 = P 2 + P 2 U * 2 , 11 

and (l/2)u* 2 (l/2)u* 2 G q 

(E2 + P2V2+ 2 ~^— ) - (Cl+PlVl+ 2 ~^~ ) + = °- 111 

P2U* 
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It is convenient to define 



q = G 0 / p2U 2 

where q is the specific energy which the element of mass has acquired 
from the laser pulse as it moves through the absorption region, see 
Figure 5. 

(U) So far only three relations have been specified relating four 
unknowns, but a number of useful results can be obtained before the 
equation of state of a specific material is introduced. The equations 
are most easily manipulated when written in terms of V. With this modi- 
fication, Eq. I becomes 

U*/U* = Vj/V 2 , 

1 2 



and Eq.s I and II can be manipulated to obtain 

i 



j 




P 2 -Pi 

vT^vT 



u*/Vi 

i 



CD 



These relations are of special use when ui = 0 (ambient air at rest); so 
that u* = -D. With this additional condition 

l 



| D h V '\E- 



( 2 ) 



this is the speed of the discontinuity in the LAB coordinate system. 

(u) It is possible to show that uj < D if the discontinuity is to 
propagate as indicated, so u* must always be negative. But Vi must be 
positive; so 






(P 2 -Pi)/(Vi-V 2 ) 



must be intrinsically negative. This factor will always be assumed to be 
negative. 

(U) When ui = 0, Eq. I can be written (use Eq. 2) 
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U2 




(P2-Pl)/ (V1-V2) 



( 3 ) 



N.B. that U 2 is in the LAB system. U 2 < D for all possible discontin- 
uities (otherwise the material behind the discontinuity would catch up 

to it!), so u* must also be negative. 

2 

(U) It is convenient to remove the microscopic variables u* 2 and 
u* 2 from Eq. III. From Eq. I 

u * 2 “ u * 2 = -(V1+V2MP2-P1), ( 4 ) 

and Eq. Ill becomes 

£2-ei + 1 / 2 (V 2 -Vi) (P1+P2) = q ( 5 ) 



The negative sign in the definition of q was chosen so that Eq. 5 could 
be interpreted physically. It is consistent with the fact the u* is 
intrinsically negative; so q as defined will always be positive. 



(U) To proceed further one must assume an equation of state for 
the material. The ideal gas will be used in this discussion because it 
is sufficient for demonstrative purposes, even though it can hardly be 
expected to describe a plasma correctly. In terms of the specific heat 
ratio (0^ is the specific heat at constant pressure, at constant 
volume) 



Y 2 C p /C v , ( 6 ) 

the ideal gas law can be written as 

P = (y-l)pe. (7) 

At times it is more convenient to define a temperature 
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T = e/C 



v 



( 8 ) 



(note that C y is assumed to remain constant); then 



P = (y-l)pC v T. 



(9) 



A useful auxiliary relation is 



( ^rr ) pv = e + PV = h = C p T 



( 10 ) 



(U) These basic relations will now be applied to various LSA wave 
propagation conditions. 

4.2 ONE-DIMENSIONAL LSD WAVES 

(U) Detonation waves occur when the speed of the wave exceeds the 
speed of sound in the ambient air, D > ci. The wave propagation is not 
stable unless it also meets the condition U 2 £ C 2 £ D, where C 2 is the 
speed of sound behind the discontinuity. With these conditions the dis- 
continuity propagates as a shock wave, and many of the analytic results 
obtained from shock wave studies also apply to LSD wave propagation. 

(u) The system is easiest to analyze if the CHAPMAN-JOUGET condition 
is applied: 



u* = c 2 . 



2 



which, when combined with the ideal gas medium, gives 



u * 2 = c 2 = Y2P2V2. 



(ID 



2 2 



Zel'dovich and Raizer 31 discuss the meaning of this condition in con- 
siderable detail. 



(u) When Eq. 11 is used with Eq. II, one obtains 
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P 2 = 



PiVi + u* 2 

Vi(Y 2 + 1 ) ; 



( 12 ) 



in conjunction with Eq. j t it gives 



V 2 = ( — 



u* 2 + PiVi y 2 



Vi 



)C 



Y2 + 1 



) , 



(13) 



and in Eq. Ill with Eq. 10 it reduces to 



Y 2 Yi , u*“ 

( ) P 2V 2 - ( ) p i V i - J V2P2V2 - = q 



(14) 



Note that the change in medium properties behind the shock front has 
been accounted for in this analysis by allowing Yi ^ Y 2 • If one now 
eliminates P 2 and V 2 from Eq. 14, the resultant equation is a bi-quadratic 
in u* 2 . The solution to this equation, with the replacement 



P 1 V 1 + (Yi - l)C vi Ti 



is 






1/2(Y2-1) [(Y2+l)q + (Yi+Y 2 )C Vi T 



1/2I 

* J 



jyi/2(Y 2 +i)[(Y2-i)q + (y 2 -yi)c Vi Ti] ^ 



(15) 



This equation becomes very simple in the strongly illuminated condition 
where q > > C y Ti (the specif id internal energy of the ambient gas is 
negligible compared to the specific energy added by the laser). It 
becomes 



ju*/ = i D j = jW - l)qj > 



(16) 
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Now restate Eq. 12 to provide the pressure ratio, 

P2 “t 2 ♦ (Vi - dc vi tx 
p 7 = (Y2+Yl)(Yl - DC^D * 

and Eq.s 13 and 17 provide the volume ratio 



( 17 ) 



V 2 

vT 



Y 2 

y 2 +l 




+ (Yi - DC Ti 

Vi 




(18) 



If one now uses the requirement (Yi - 1)C Ti 



< < u* 2 , Eq. 10 reduces to 



V 2 Y 2 
VT " Y2+1 



(19) 



which applies in the strong shock limit of a strongly fed wave which also 
satisifies the CHAPMAN- JOUGET condition. To the same approximation 



P 2 2 (Y 2 -I) q _ YxD 2 

PT = (Yi-D C Vi Tj = (Y2+1)C 2 

P 2 Yi 

* Y 2 +I M C20) 

where M is the Mach Number of the detonation wave motion. 

(U) If one replaces q by its definition in Eq. 16, one obtains a 
relation between the wave speed D and the irradiance of the laser: 

1/3 

D = [2(y* - l)ViG 0 ] (21) 
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Note that this result depends upon the ambient density in front of the 
shock, the specific heat ratio behind the shock, and the irradiance. 

The relation between D and G has been well tested experimentally. 

4.3 ONE -DIMENSIONAL LSC WAVE 

CuD Combustion waves occur when the velocity of propagation of the 
disturbance is less than the velocity of sound. Velocities for LSC 
phenomena range from D = 0, the plasmatron, to the order of a few 
meters/second. Here the D = 0 case is specifically excluded. 

(U) The LSC disturbance is a more complex structure than is the 
LSD wave as sketched in Figure 6. This figure, from the work of Jackson 
and Neilsen 32 , is the result of a computer computation. Its qualitative 
features agree very well with the experimental data available on LSC and 
LSD waves. 

(U) The LSC wave structure is conveniently divided into three 
regions: A, the ambient air in front of the disturbance; B, the transi- 

tion region, which has been warmed and which contains some ionized par- 
ticles, but which is transparent to the laser photons; and C, the plasma , 
which is opaque to the laser photons. The electron density and tempera- 
ture increase continuously with distance from the wave front so that at 
the plasma front the plasma frequency of the ionized medium equals the 
frequency of the incident laser light. 

(U) The LSC wave differs from the LSD wave because the disturbance 
created by the wave front is not sufficiently intense to produce the 
electron density required to support the plasma; so additional ionization 
must be generated by the plasma. For the LSD wave the two fronts coalesce 
and the analysis is much simpler. 

(U) Because of the more complex structure, the simple conservation 
law analysis of the LSD system contained in Eq.s I } 11^ and III is no 
longer appropriate. It could be applied between regions A and C, but 
most of the interesting phenomena which control the wave development 
occur in region B. In region B the parameters are functions of distance 
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TEMP 

-12 

(10 etgs ) 



PRESSURE 

(atm) 



TEMP 

(lO 1 ^ etgs) 




Fig. 6A One -dimensional LSD wave driven by a CO^ laser flux 

of 1 x 10 7 watts/cm 2 decoupling of high temperature region 
and shock wave does not occur = 4.0 x 10 s cm/sec. 
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Fie. 6B One-dimensional LSC wave and shock wave precursor 
driven by a C0 2 laser flux of 1 x 10 6 watts/cm 2 . 

V LSC = 1,1 x 1C>5 cm / sec an<J V SH0CK = 1,4 x 105 cm / sec 
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from the wave front; so a full divergence theory analysis is required. 

In addition, the functional variations are driven by the atomic level 
process discussed in Section 2; so the resultant non-linear differential 
equations are extremely complicated. 

(u) Here it seems appropriate to discuss the model and the pheno- 
menological parameters, but a presentation of the actual numerical 
methods used to solve the equations would go beyond the intent of this 
survey. The interested reader will be referred to the current literature 
at the appropriate points in the discussion. In many respects this 
presentation draws heavily upon Nielsen and Canavan 5 , but there are 
significant variations from their notation. The original analysis of 
LSC waves performed by Raizer 33 , has strongly influenced all subsequent 
studies. 

(U> Because the LSC wave propagates so slowly, to a very good 
approximation the pressure may be assumed constant throughout the dis- 
turbance and the conservation of momentum condition contributes nothing 
to the analysis. Similarly, there is very little mass flow; so effec- 
tively only the energy conservation condition remains. This condition, 
written in differential form is 

Divfpvh + Pv) = Q. 

Here h is the enthalpy and Q contains ail of the energy source and sink 
terms. Three terms can be identified in Q: 

Q = $ + Y + E, (22) 

where <J> is the radiation term, T the heat conduction term, and £ is the 
electron transport term. Each of these terms operates somewhat differently 
in the three regions. Q vanishes in region A, because by definition the 
medium is undisturbed there. The laser supplies energy in region C, 
because the plasma is opaque, so $ is a source term which can be written 

* = K o G o (x) (23) 

where K q is the absorption coefficient, and G(x) is the residual laser 
flux at distance x -from the wave front which is taken as the origin. The 
laser intensity in turn must satisfy 
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dG 

dx 



= -KG 



(24) 



subject to the boundary condition G(o) = G the laser irradiance. 

o j 



(U) The absorption coefficient, K, is a strongly varying function 
of temperature, see Figure 7. Because the coefficient rises so sharply, 
in most calculations it is approximated by the function 



K(T) =0.0 T < 12, 000°K 

= 0.7 cm" 1 . T > 12,000°K (25) 

Other investigations prefer to set K(T) = 0.6 cm" 1 when T > 1.0 eV 32 . 



(u) In region C, the thermal conduction term, ¥, is a loss term. 

It supplies much of the energy which increases the temperature of 
region B. The thermal conductivity of air, 8, is an extremely compli- 
cated function of temperature, see Figure 8. A convenient approximation 
for 8 is 18 

8 = 2 x 10~ 6 T (watts/cm °K) . (26) 



The heat flux is given by 

d 2 T 

Y = -8 — . 
dx 2 



(27) 



To a reasonable approximation the temperature of the plasma is constant, 
Tp ~ 1 eV; so there is heat conduction from the plasma into region B. 

By analogy, with chemical combustion waves, Raizer 33 felt that heat 
conduction would supply sufficient energy to increase the temperature of 
region B to "combustion". Careful calculation by Hall, et.al. 18 and 
Thomas, et.al. 28 have shown that this is not the case. 



(U) Jackson and Nielsen 32 have investigated the heating of region 
B by radiation transport from the plasma of region C. The opaque 1 eV 
temperature plasma emits strongly throughout the visible region of the 
spectrum and also well into the ultraviolet. The dominant radiation 
heating contribution in region B comes from this source rather than 
from the laser beam. Effectively, the plasma serves as a frequency 
converter for the laser energy, and re-emits it in a spectral region 
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ABSORPTION COEFFICIENT (CM 




*-T(K) 



Fig. 7 Abortion coefficient of a carbon dioxide 
laser in air as a function of temperature. 

From reference 18. 
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Fie. 8 Thermal conductivity of air as a function of 
temperature. From reference 18. 
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where the weaker plasma of region B is already opaque. 
(U) The radiant flux term is written as 



$ = -c \ dvK* (G (T) - G ) (28) 

J vp v J K J 

00 

in the Jackson and Nielsen model. Here G^ is the irradiance per unit 
frequency from all sources and G (T) is the equilibrium radiative 
energy flux density. In this formulation $ is a loss term in region C 
and a source term in region B. The simple attenuation approach of 
Eq. 23 comprises one part of this expression. The coefficient is the 
absorption coefficient at frequency v corrected for induced emission, 
see Zel'dovich and Raizer 31 . 

(U) This coupling of the various portions of the LSC wave through 
radiation transport makes the entire model strongly non-linear. Still, 
the model equations can be solved analytically if three assumptions are 
made : 

1) The absorption coefficient must be treated as a step function 

as in Eq. 25. (Jackson and Nielsen use 0.6 cm 1 with a threshold tempera- 
ture of 1 eV) . 

2) The ratio B(T)/Cp(T) is set equal to a constant f = 1.0 x 10” 3 
gm/cm-sec. The enthalpy of air is shown in Figure 9, together with a 
useful analytic approximation. The relation = (Bh/BT^) and this 
approximation give f ^ 2 x 10 3 gm/cm-sec. 

3) The radiative transport term $, is written as an exponential 
heating term in front of the plasma front and a constant absorption term 
behind the plasma front. 

(u) The system of equations can then be solved by a series of 
iterations which adjust the values of the necessary parameters. Jackson 
and Nielsen find that the most important contribution energy contributes 
from radiation transport comes in the range 14.5 eV <_ V <_ 20 eV, which is 
well into the uv region. Account must be taken of the anisotrophy of the 
radiation field. The authors use 19 frequency groups to establish the 
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TEMPERATURE 




Fig. 9. Enthalpy of air as a function of temperature. 
From reference 18. 
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self-consistent solutions. 

(u) The model preduces the correct trend of the propagation 
velocity, D, with the laser irradiance, 3, but agreement in magnitude 
is not satisfactory; see Figure 10. 

(U) Jackson 34 has shown that the presence of the target surface 
forces a streaming velocity upon the air in which the LSC wave is propa- 
gating; so that the experimentally determined velocity is considerably 
larger than the radiation transport model predicts. 

(u) The final transport term, £, the electron fluence, has appar- 
ently not been considered by any investigator. Clearly region C contains 
a high density of electrons and no free electrons are present in region 
A; so there must be a diffusion of electrons into region B driven by a 
concentration gradient. The effect is complicated by the fact that ions 
must accompany the electrons to maintain charge neutrality; so the model 
is one in which a high density plasma diffuses into a low density plasma. 
This type of an analysis would require an additional conservation law, 
the number of particles, and a kinetic theory approach to the spatial 
distribution of particles. 

(U) One other important idea is suggested by Figure 6, which 
depicts both LSC and LSD waves produced by the same computer code. With 
the LSC wave irradiation condition, the LSC wave, the plasma front, falls 
behind the shock front; so the wave presents the characteristic double 
structure. As the irradiance increases, the speed of the plasma front 
increases until finally the two front propagate together as an LSD wave. 
Apparently there is a continuous variation from plasmatron through LSC 
behavior to LSD behavior. 

(U) Jackson 35 has calculated the irradiance at which the LSD behavior 
should appear. He finds 



G t = [ 



(y-i ) 1 

2y 



J p A h B 



3 2 



where is the ambient air density and h fi is the enthalpy in region B. 
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KLOSTERMAN AND BYRON DATA 
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Fig. 10 Plasma wave velocity as a function 
of irradiance. 
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With reasonable values for the parameters, H = 2 x 10 7 W/cm 2 which agrees 
quite well with the experimental results of 2.25 x 10 7 W/cm 2 18 . 

4.4 LSA WAVES IN TWO-DIMENSIONS 

(U) The geometry of any LSA wave system is clearly cylindrical as 
soon as it enters the transition phase. The one-dimensional analyses 
presented here are only approximations which have meaning when radial 
effects (across the wave) can be neglected relative to longitudinal 
effects (parallel to the direction of propagation). When ignition, trans- 
ition, and propagation are under study radial effects introduce additional 
loss mechanisms which withdraw energy and electrons from the plasma. 

(U) Magnetic confinement of the plasma, caused by its motion, will 
be neglected in this discussion. It is an additional consideration which 
will tend to reduce the effect of the radial loss mechanisms. 

(u) The significance of the radial direction is effectively deter- 
mined by the diameter of the laser beam in the volume under study. A 
small diameter beam allows significant radial transport of electrons, 
heat, and radiation; a large diameter beam does not. Raizer 3 has shown 
that when the amount of radial of expansion is small the energy balance 
equation should be corrected by the replacement G -* G/ (1 + £/R) , where 
R is the radius of the laser beam, and L is the length of the absorption 
region; that is the depth of the plasma in region C of the last section. 
Jackson and Nielsen 32 have done rather careful calculations which confirm 
that radial losses are of secondary importance in the analysis of LSA 
waves when the condition t < < R is satisfied. 
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5.0 SUMMARY 



(U) This article reviews the initiation and propagation of LSC and 
LSD waves. It concludes that prompt ignition requires that electrons be 
emitted from the target surface, most probably at physical features 
and/or inclusions in the surface. The electrons are followed by ions 
extracted from the surface by electrostatic forces. If the irradiation 
is sufficiently intense the plasma moves away from the surface and feeds 
upon the incident laser radiation. The plasma is accompanied by a 
shock wave. At intermediate irradiance levels the plasma falls behind 
the shock wave and an LSC wave results. For irradiance levels above 
~ 2 x 10 7 W/cm 2 the plasma and shock wave move together as an LSD wave. 
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